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Gas-phase multiply charged proteins have been sympathetically cooled to translational tempera- 
tures below 1 K by Coulomb interaction with laser-cooled barium ions in a linear ion trap. In one 
case, an ensemble of 53 cytochrome c molecules (mass :^ 12390 amu, charge +17 e) was cooled by 
^160 laser-cooled barium ions to less than 0.75 K. Storage times of more than 20 minutes have 
been observed and could easily be extended to more than an hour. The technique is applicable to 
a wide variety of complex molecules. 



In the last years, the development of techniques for 
producing cold molecules (< IK) has moved into 
focus. Remarkable advances in this young field of 
quantum physics are enabling studies of light-molecule, 
atom-molecule and molecule-molecule interactions in a 
temperature regime not accessible previously. Cold 
molecule types produced so far range from neutral di- 
atomic molecules produced by photoassociation in ultra- 
cold atomic gases [1] to polar few-atom molecules slowed 
down by Stark deceleration [2] and charged diatomic to 
polyatomic molecules sympatheticallv cooled by laser- 
cooled atomic ions in ion traps (3, H HI- The heaviest 
cold molecular species prepared had a mass of 410 amu. 

It is of interest to develop methods that allow a similar 
control over much heavier particles like anorganic clus- 
ters or biomolecules, consisting of hundreds or thousands 
of atoms. So far, such complex systems have been stud- 
ied in ion traps using buffer gas cooling [g*, 7] . The buffer 
gas thermalizes with a cryogenically cooled chamber con- 
taining the trap and collisionally cools the trapped ions 
internally and translationally to temperatures of ~ 10 K 
[S]. Recently we have reported on a technique for cool- 
ing polyatomic molecular ions that allows one to achieve 
significantly lower translational temperatures and long 
storage times of up to hours. As examples, singly pro- 
tonated dye (mass 410 amu) Q and glycyrrethinic acid 
molecules (mass 471 amu) [9] produced by electrospray 
ionization (ESI) have been sympathetically cooled to less 
than 150 mK by their mutual Coulomb interaction with 
laser-cooled ^^^Ba^ ions in a linear quadrupole trap. 

Here we demonstrate a major extension of this tech- 
nique to much heavier molecules (such as proteins), tak- 
ing full advantage of the ability of the ESI ion source to 
produce multiply charged ions with mass-to-charge ratios 
suitable for conventional quadrupole mass filters and lin- 
ear ion traps jiol, [ll|. We have produced, trapped, sym- 
pathetically cooled, and maintained cold for extended in- 
tervals 12- and 17- fold protonated protein molecules (cy- 
tochrome c, a protein consisting of about 100 amino acids 
with a total mass of 12390 amu) using Doppler laser- 
cooled ^^^Ba+ ions in a linear quadrupole trap. The two 
trapped species were detected and distinguished via ex- 
citation of their specific motional resonances. The sym- 



pathetic cooling was qualitatively proven by measuring 
the energy distribution of the ions and quantitatively an- 
alyzed using molecular dynamics simulations. 

Our apparatus ^ is shown schematically in Fig. [TJ It 
consists of an ESI ion source for the production of singly 
and multiply protonated molecular ions, a quadrupole 
mass filter for the selection of specific molecular species, 
a rf octopole ion guide to transfer the selected molec- 
ular ions from the medium vacuum region of the ap- 
paratus to an ultrahigh vacuum (UHV) chamber, and 
a linear quadrupole trap in this UHV chamber to store 
the molecular ions for cooling and further investigations. 
The preparation of the laser-cooled ^^^Ba+ ion ensembles 
used for sympathetic cooling of the molecular ions as well 
as the required lasers (493 and 650 nm) is described in [1] . 
With sufficient cooling power the barium ions arrange in 
ordered structures, so-called Coulomb crystals, that can 
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FIG. 1: (Color online) Schematic overview of our setup. 
Singly or multiply protonated molecular ions from an ESI 
ion source are selected by a quadrupole mass filter and trans- 
ferred via a RF octopole ion guide to a linear quadrupole trap 
in an UHV chamber, (a) CCD image of a typical trapped 
barium ion Coulomb crystal used as coolant for the molecu- 
lar ions, (b) Secondary structure and ESI mass spectrum of 
cytochrome c showing the distribution of the different grades 
of protonation. 
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be imaged with a charge-coupled-device (CCD) camera 
[Fig. [1] (a)]. The images provide information not only 
on the fluorescing ^^^Ba+ ions, but indirectly also on the 
numbers and temperatures of the invisible sympatheti- 
cally cooled molecular ions. Molecular dynamics (MD) 
simulations are used to derive these data from structural 
deformations and sharpness of the barium ion Coulomb 
crystals [12]. 

Multiply protonated molecules of cytochrome c from 
horse heart (Fluka BioChemika), (cyt c + nH)^+ with 
grades of protonation, n = 7 — 17, have been produced 
by ESI from a solution concentration of 10~^M in 1:1 
methanoliwater. As the maximum n depends on the pH 
value of the solution, a solution with a higher concentra- 
tion of acetic acid of 2% was used to produce 17- fold pro- 
tonated cytochrome c molecules, here denoted as Cyt^^+, 
and a lower concentration of 0.5% was used to obtain 12- 
fold protonated molecules Cyt^^+. 

The preparation procedures for sympathetically cooled 
Cyt^^+ and Cyt^^+ ion ensembles are identical. First, 
helium buffer gas at room temperature is injected into 
the trap chamber to pressures of 1 — 10 x 10~^ mbar in 
order to subsequently reduce the kinetic energy of the 
molecular ions from the octopole ion guide so far that 
they can be trapped. After loading ions of the desired 
molecular species for 1 — 10 s depending on the desired 
ion number and the ion flux (typically several thousand 
ions per second reach the trap), the buffer gas is pumped 
away and the background pressure in the trap chamber 
returns to < 1 x 10~^ mbar within a few minutes. Then, 
barium ions from an evaporator oven are loaded into the 
trap and cooled by the two cooling laser beams propa- 
gating along the trap axis. Impurity ions generated dur- 
ing loading (such as molecular fragments or CO2" from 
the evaporator) can be removed from the trap by ap- 
plying an additional appropriately strong ac voltage to 
the trap electrodes in order to excite the ions' specific 
motional resonances beyond stability [l^. At a typical 
rf amplitude of 450 V the radial resonance frequency of 
^^^Ba+ ions is 118 kHz, while that of cytochrome ions is 
in the range of 10—22 kHz depending on their protona- 
tion grade. At the end of the preparation procedure, two 
frequency scans (125—1000 and 110—40 kHz) are applied 
to purify the barium-cytochrome ion ensemble. 

The same principle of radial excitation can be used in 
a nondestructive way to identify trapped ion species. For 
this purpose, frequency scans with ~ 10 times lower am- 
plitudes than for purification are performed while acquir- 
ing the barium ions' fluorescence intensity with a photo- 
multiplier. Under resonance excitation of a particular 
ion species, its motion heats all trapped ions due to their 
mutual Coulomb interaction. As a consequence, the bar- 
ium ions' fluorescence drops, leading to dips in the ra- 
dial excitation spectrum at the species-specific resonance 
frequencies [15I. Il6l|. Cyt^^+ and Cyt^^+ in barium- ion 
Coulomb crystals have been identified and distinguished 
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FIG. 2: Radial excitation spectra of barium crystals con- 
taining differently protonated cytochrome ions. The upper 
(Cyt^^^) and lower (Cyt^^^) curves were acquired at the 
same trap parameters and yield the same barium frequency 
of 95.4 kHz used for calibration. The obtained cytochrome 
frequencies are well confirmed by their theoretical values. 

via their radial resonance frequencies as shown in Fig. [21 
The upper and lower curves are composed of two differ- 
ent scans: one from 5 to 65 kHz for the detection of 
the molecular ions and one from 40 to 150 kHz with 
a 10 times lower amplitude for the barium ions (in or- 
der to avoid crystal destruction). As the rf amplitude 
was the same for all four scans (362 V), the barium 
resonances coincide at z/^(Ba+) = 95.4 kHz and can 
be used for calibration. With the resonance frequency 
^ Q/tti (for —a <C q [13]) the expected frequen- 
cies for (noninteracting) cytochrome ions can be calcu- 
lated as i/^(Cyt^+) = i//(Ba+) • [Q(Cyt^+)/Q(Ba+)] • 
[m(Ba+)/m(Cyt^+)], which yields 12.8 kHz for Cyti2+ 
and 18.1 kHz for Cyt^^+. These values are well con- 
firmed by the experimental results of 13.5 and 18.6 kHz. 
The small differences of less than 5% can be ascribed to 
various factors like the direction and amplitude of the 
frequency scans [15] or motional coupling effects [16]. 

Additionally, we have implemented a destructive mass 
identification technique which is based on a mass- 
selective extraction of the ions from the trap. When de- 
creasing the trap rf amplitude t/rf the trapping stability 
parameter q [13] is reduced and ions escape the trap at a 
mass-to-charge ratio dependent amplitude U^l ex ^Jm|Q 
. In our setup an ion detector below the trap registers 
the escaping ions during such a controlled reduction of 
the trap rf amplitude, leading to an ion extraction mass 
spectrum as shown in Fig. [3l Although this technique is 
destructive and less accurate than the radial excitation, 
it has several advantages. It is fast (an extraction takes a 
few seconds) and it works also with uncooled ions. More- 
over, it provides qualitative information on the ion tem- 
peratures as the mass spectra of warmer ion ensembles 
are broadened due to the ions' wider energy distribution. 
This is shown in Fig. [3] for Ba+/Cyt^^+ ensembles in 
the cooled (black) and the noncooled case (gray). In the 
cooled case, both the peak of the laser-cooled Ba+ ions 
and the peak of the Cyt^^+ ions are narrower. This is a 
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FIG. 3: Ion extraction mass spectra of Ba^/Cyt^^^ ensembles 
at different temperatures. In the noncooled case (gray) the 
mass spectrum is broadened compared to the laser-cooled case 
(black). Here, not only the peak of the laser-cooled Ba^ ions, 
but also that of the simultaneously trapped Cyt^^^ ions is 
narrower, which is a proof of their sympathetic cooling. 

qualitative proof of their sympathetic cooling. 

Ba+/Cyt^+ crystals were routinely produced. Several 
Ba+/Cyt"'^^+ and Ba+/Cyt"'^^+ ion crystals have been an- 
alyzed quantitatively by the use of MD simulations in or- 
der to obtain the numbers of the trapped ions and their 
translational temperatures Tsec,i = where {Ei) 

is the time- and subensemble-averaged secular kinetic en- 
ergyper ion of the species i. Our MD method is described 
in Q and Here, we do not consider micromotion 

in the simulations explicitly since it does not affect the 
CCD images [l^. Micromotion-related effects, such as rf 
heating, are taken into account phenomenologically. 

In a quadrupole trap ions of different charge-to-mass 
ratios arrange radially separated according to their ef- 
fective radial trap potentials <l>(r) ex [13i] and the 
interspecies repulsion ^ Q1Q2' In the limit of cylindri- 
cal symmetry, this radial separation scales as ri/r2 = 
((52^1 /Qi^2)^^^ with the outer radius ri of the lower 
mass-to-charge ratio (mi/Qi) subensemble and the in- 
ner radius r2 of the higher mass-to-charge ratio (m2/Q2) 
subensemble [l9^. As shown in Figs. [4] (a) and[4](d), the 
multiply charged cytochrome ions form a sheath around 
the barium ions which radially squeezes and axially pro- 
longs the barium ion subensemble. In the experiment two 
CCD images of the ion crystals are recorded: one with the 
molecular ions present and one after the molecular ions 
have been removed. In the case presented in Figs. [41(a) — 
[4] (c) the Cyt^^+ ions prolong the barium ion subensem- 
ble axially by 240 jam [Fig. [41(b) and [4](c)], which can be 
explained by the presence of 53 Cyt^^+ ions, according 
to MD simulations. The barium ensemble consists of 158 
laser-cooled ^^^Ba+ ions and 230 barium isotopes other 
than -'^^^Ba. Due to nonoptimal trap parameters during 
the experiment of Fig. [4] (b), the observed radial squeez- 
ing of the barium ion subensemble is less pronounced 
than in the simulated image. 

For the temperature determination a specific heating 
rate hi = d{Ei)/dt for each ion species i is chosen that 
phenomenologically describes heating effects such as col- 





FIG. 4: (color online) Experimental and simulated CCD im- 
ages of a Ba^/Cyt^^^ and a Ba^/ Cyt^^^ ion crystal, (a, d) 
Simulated images in radial (left) and axial view (right), as 
they would appear if all ions would fluoresce. The cooling 
lasers propagate to the left, separating the ^^^Ba^ ions (blue) 
and the barium isotopes (red) axially due to the light pressure 
force. The cytochrome ions (green) form a sheath around the 
barium ion subensemble. (b, e) Experimental image of the 
barium-cytochrome ion crystal and its simulation, showing 
only the fluorescing ^^^Ba^ ions, (c, f) Crystals after the re- 
moval of the cytochrome ions. Here, the barium ion ensemble 
is no longer deformed nor heated by the molecular ions. 



lisions with background gas particles, rf heating, or elec- 
tric field noise [12]. For the laser-cooled ^^^Ba+ ions 
a cooling rate is employed corresponding to the laser 
cooling force Fd = —Pvz^ where P is the damping co- 
efficient and Vz the z component of the ion velocity 
with the cooling lasers propagating in the positive z di- 
rection [5"]. The other species are cooled sympatheti- 
cally through the interparticle Coulomb interaction, re- 
sulting in a steady state in which every ion species i 
has its specific translational temperature Tsec,i- For the 
pure Ba+-ion crystal shown in Fig. [H (c) a heating rate 
^Ba+ = f /cb(3.7 K/s) [2^ needs to be applied to repro- 
duce the experimental CCD image concerning shape and 
internal structure of the ^^^Ba+ ion subensemble. Here, a 
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realistic damping coefficient /? = 1.7xl0 ^^kg/s was as- 
sumed [21], which yields the temperatures Tg^^ i38Ba+ = 
19 mK and T^^c ,Ba+ isotopes ~ mK. As all Condi- 
tions remained constant, the same damping coefficient 
(3 and the same heating rate h^^+ are applied to re- 
produce the CCD image of the corresponding "parent" 
Ba+/Cyt^^+ ion crystal [Fig.[4](b)]. Here, a heating rate 
hQy^i7+ = |/cb(71.6 K/s) must be present in order to 
reproduce the observed warmer barium ion subensem- 
ble whose temperatures are Tg^^ i38Ba+ = 
^sec,Ba+isotopes = 109 mK. With this value, the cy- 
tochrome temperature is T^ec,Cyt^'^+ — '^04 mK. As the 
laser cooling efficiency varies depending on numerous fac- 
tors, we consider a range of (5 values, as observed by 
[22! . To obtain a lower (upper) limit for the heating rates 
and temperatures we simulate the cases = 1.0(2.0) x 
10~^^ kg/s. For the ion crystal shown in Fig. [4] (b) with 
^sec,i38Ba+ = 66 mK, this yields the following values: 
W = f ^B(3.7i?:^ K/s), Tg,,,Ba+ isotopes = 1091?^ mK, 
hcyt^i^ = 1^5(71.6129.2 K/s), and Tg^c,Cyti^+ = 
704t^^ mK. 

The temperature of the molecular ions mainly depends 
on the laser cooling power, the pressure in the vacuum 
chamber, and the ratio of the numbers of laser-cooled and 
sympathetically cooled ions and also their mass-to-charge 
ratio compared to that of the laser-cooled ^^^Ba^ ions. 
The closer the ratios are, the lower is the radial sepa- 
ration of the two species in the trap and the higher is 
the sympathetic cooling efficiency. In the cases shown in 
Fig. m the radial gap (r2 — ri) between the Cyt^^+ and 
the ^^^Ba+ ion ensembles is 64 //m and in the Cyt^^+ case 
87 /im. This contributes to the fact that the Cyt^^+ ions 
could only be cooled to 1.2 K, although the barium ions 
at 26 mK were much colder than in the Cyt^^+ case 
[2^ 1. As a consequence of the different cytochrome tem- 
peratures, the Cyt^^+ subensemble [Fig. [H (d)] appears 
more blurred in the simulated images than the colder 
Cyt^''+ subensemble [Fig. [J] (a)]. 

In conclusion, we have shown that heavy and thus com- 
plex particles, appropriately charged, can be translation- 
ally cooled in the gas phase to temperatures significantly 
below those reachable with other techniques. The lowest 
temperature observed was below 0.75 K. The tempera- 
ture is limited by heating experienced by the molecular 
ions. Its origin is uncertain so far, but may be due to 
electrical noise. Combining the present result with ear- 
lier work (sympathetic cooling of masses up to 410 amu) 
[1, m we reasonably expect that ions of any mass be- 
tween 32 and 13000 amu can be sympathetically cooled, 
provided they have an adequately high charge. The par- 
ticular advantage of this cooling method is that the ions 
can be stored for more than 20 min (easily extendable to 
hours using an improved vacuum to avoid loss reactions 
of the ^^^Ba+ ions) with a strong spatial localization in a 
well-defined and nearly collisionless environment. These 



are ideal conditions for, e. g., studies of low-energy pro- 
cesses, like cold collisions or triplet-singlet decay, or high- 
resolution spectroscopy for structural or conformational 
characterizations of biomolecules. For many experiments 
cooling of the internal degrees of freedom of the molecules 
would be favorable and could be implemented by radia- 
tive cooling in a cryogenic environment. 
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